Slow relaxation in microcanonical warming of a Ising lattice 
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Abstract. Wc study the warming process of a semi-infinite cylindrical Ising lattice initially ordered and 
coupled at the boundary to a heat reservoir. The adoption of a proper microcanonical dynamics allows a 
detailed study of the time evolution of the system. As expected, thermal propagation displays a diffusive 
character and the spatial correlations decay exponentially in the direction orthogonal to the heat flow. 
However, we show that the approach to equilibrium presents an unexpected slow behavior. In particular, 
when the thermostat is at infinite temperature, correlations decay to their asymptotic values by a power 
law. This can be rephrased in terms of a correlation length vanishing logarithmically with time. At 
finite temperature, the approach to equilibrium is also a power law, but the exponents depend on the 
temperature in a non-trivial way. This complex behavior could be explained in terms of two dynamical 
regimes characterizing finite and infinite temperatures, respectively. When finite sizes are considered, we 
evidence the emergence of a much more rapid equilibration, and this confirms that the microcanonical 
dynamics can be successfully applied on finite structures. Indeed, the slowness exhibited by correlations 
in approaching the asymptotic values are expected to be related to the presence of an unsteady heat flow 
in an infinite system. 



1. Introduction 

The transient regime of an Ising system starting from equilibrium at a uniform temperature and then 
suddenly coupled to a thermostat at a different temperature, is an issue involving both practical problems 
and theoretical aspects of non equilibrium. For instance, as it is well known, cooling a system at a 
temperature smaller [I] [3J El H] or equal [3 E] to the critical one, gives rise to interesting dynamical 
out-equilibrium phenomena, where the transient regimes become indefinitely long. This approach, even for 
an artificial Monte Carlo dynamics, provides a deep insight in important physical problems, such as aging 
phenomena and coarsening of magnetic domains. 

We remark that within the usual canonical approach used in these studies (e.g. Metropolis), the 
interaction with the reservoir yields a uniform thermalization, a fact avoiding heat transport. Therefore, in 
a canonical frame, the investigation of heat propagation during a out-equilibrium warming process is not 
even possible. This excludes interesting issues related to transport: which are the relevant parameters of 
the warming process, the role of boundary conditions, the space-time scales, the evolution of correlations 
aligned or transversal to the heat flow, etc. [7J Ul |H1 HP] . In this canonical context heat transport can be 
studied only at interface separating two systems at different fixed temperatures [llj . 

Thus, in order to investigate heat transport phenomena in a transient regime, the appropriate 
dynamical approach should be the microcanonical one. This approach allows to couple the system to 
a thermostat and to investigate how the conserved energy flows within the system itself. However, so far, 
this approach has been largely neglected, possibly because of the need of a unique dynamical frame in 
a wide range of temperatures and/or of underlying geometries, not matched by standard microcanonical 

rules [una muslin nam]. 

Now, a new efficient microcanonical dynamics has been recently introduced in [19] . allowing for 
enhancements in the study of heat flow for a spin lattice coupled to thermostats at arbitrary fixed 
temperatures. Indeed, among the advantages of this dynamics (including e.g. the simple definition of 
local temperatures, the direct computability of conductivity, the possibility of applying to a great variety 
of discrete topologies and coupling patterns), there is also a property that the previous rules lacked, i.e. 
an efficient behavior in the whole range of temperatures from to oo. We remark, that previous results 
for such a dynamics on homogeneous or disordered systems, are limited to the stationary regime, which 
is established after a long-time interaction between the systems and the thermostats [TH1 HOI 121] ■ On the 
other hand, as already mentioned, here we focus on genuine non-stationary phenomena emerging when the 
system is suddenly put in contact with a thermostat at different temperature. 
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In particular, the system under study is an initially ordered rectangular lattice, where at time one 
side is coupled to a reservoir at fixed temperature T, and consequently a heat diffusion takes place in the 
direction orthogonal to the hot side. This rectangle, in order to have a permanent state of transience, 
should be thought of as infinite in the diffusing direction. In practice, we shall take this size L greater 
than the grid thermalized by the expanding heat front during a long observation time. Periodic boundary 
conditions are assumed in the transverse direction, so that the geometry reduces to a semi-infinite cylinder. 

Some basic concepts of out-equilibrium physics, which have been developed for instance in the study 
of coarsening process, will be very useful also in our approach. In particular, dynamical scaling is a general 
tool that on one hand gives a simple description of the system evolution, and on the other it provides 
a theoretical basis for the introduction of universal concepts, such as correlation length and dynamical 
exponents. In this perspective, since for our dynamics a diffusive behavior has been proven [THl [201 HI] , a 
generic observable O is expected to depend on the time t and on the distance x from the reservoir according 
to the dynamical scaling: 

0{x,t) ~ /(f) , t = t/x 2 . 

Notice that the rescaling procedure allows to introduce a stationary approach even in the transient problem; 
indeed, data obtained at different times can be compared to each other by a simple redefinition of distances: 

x = x/t 1/2 cx x/£(t) , 

where £(t) is a characteristic diffusive length. 

The most interesting feature we will evidence is that, in the presence of non-stationary heat 
propagation, the approach to equilibrium may be very slow, even when heating a system at high 
temperature. 

In this context an interesting observable is the correlation (cr Xty (t) o~ x> y-\- n (t)}, measured at time t 
between the value of spins at sites y and y+n of a line orthogonal to the heat flow at distance x from the hot 
side. For n = 1 this correlation corresponds to the magnetic energy. Since we are interested in the relaxation 
towards the equilibrium, we shall consider the quantity C n (x, t) = (a x ,y(t) cr X)2/ +„(i)) — (o- x ,yO~ x ,y+n) ■> where 
(o~x,yO~ x , y +n) is the correlation function of a system at equilibrium at the temperature of the heat reservoir. 
Thus, C n (x, t) tends to in the limit t — > oo; moreover, according to the scaling ansatz C n (t, x) is a function 
F n (t). 

By numerical simulations, first we check the scaling hypothesis, then we study the behavior of the 
scaling functions F n (t). In particular, after extensive Monte Carlo simulations, we obtain F n ~ t~ a ( n ' T \ 
where a(n, T) = a\{T)n + ao(T). In this context two simple behaviors emerge, namely a(n, T) ~ n/2 and 
a(n, T) 1/2, for large and small temperatures respectively. We will show that a dynamics characterized 
by a(n,T) = n/2 provides a reasonable description of the system evolution at infinite temperature. On 
the other hand, a(n, T) = 1/2 is expected to be a good behavior for finite temperature systems. In this 
perspective the complex dependence on T displayed by the exponents a(n, T) emerges as due to a dynamical 
crossover between the dynamics pertaining to finite and infinite temperature respectively. 

We notice that the diffusive behavior characterizing the dynamical scaling, as well as the exponential 
decay of the correlations as a function of n, are expected features of the model; indeed, the correlation 
decay is typically exponential in the Ising model as soon as the magnetic domains are regular, i.e. outside 
the critical regime. On the other hand, the way the correlations tend to their equilibrium values, i.e. by 
a power law in t, is not so obvious. It seems that the system persists in a transient regime even for times 
much larger than the characteristic timescale of diffusion, which is proportional to x 2 . In some way, even 
at fixed distance from the thermostat there is not a characteristic time for which the equilibrium could be 
considered as reached. 

We will finally show that such a slow relaxation is a typical feature of open infinite systems where non- 
steady heat transport is present for all times. Indeed, in a finite system, a rapid decay to the equilibrium 
value takes place as soon as the diffusion characteristic length £(t) is of the order of magnitude of the 
system size L. For this reason the dynamics considered here can be adopted for the study of the stationary 
properties of finite systems, as it was done in [Ti?l |2T)1 [2Tj . 

The paper is organized as follows: In Sec.[2]we briefly review the microcanonical dynamics exploited to 
simulate the evolution of the system, then, in Sec. [3] we present the numerical results focusing on diffusion 
features, on correlations and on finite-size effects; finally, Sec. 0] is left for conclusions. 

2. Microcanonical Dynamics and Model 

In order to be self-contained, we briefly recall the dynamics introduced in |19) . whose main novelty consists 
in assigning to the links, beyond the usual magnetic energy E™, a sort of "kinetic" term E^ which is 
a non-negative definite quantity able to compensate the variations of magnetic energy due to the spin 
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flips involving the adjacent nodes i and j. This makes an essential difference with respect to the Creutz 
microcanonical procedure jTSJ [TT] , where a bounded amount of kinetic energy is assigned to each node: in 
that case, indeed, both the energy boundedness and the connectivity of the nodes entail many limitations, 
ranging from geometrical or topological constraints to non-ergodicity at low energy density. 

Such limitations, as shown in T5] , are efficiently solved by the new dynamics. An elementary move 
for it, starting from a random distribution of link energies, consists in the following: 

(i) select randomly one of the four spin configurations for a randomly extracted link {i,j}, and evaluate 
the magnetic energy variation AE 7 ^ due to this choice. If the coupling = 1 , as in the present case, 
AE™ is an integer whose value depends on the connectivity of the nodes i and j; 

(ii) if AE"j < 0, accept the choice and increase the link kinetic energy E^ of AE^; if AE^ > 0, accept 
the choice and decrease E^ of AE™ only if the link energy remains non negative, otherwise reject the 
choice. 

The unit time step consists of N moves, where TV is the total number of links. Since a link is defined 
only by the adjacent nodes, this rule applies on arbitrary non-oriented connected graphs, providing a great 
generality to the dynamics. 

We remark that at low energies the dynamics is rather time-consuming as for most of the extracted 
links the move is typically rejected; yet, differently from the Creutz procedure, here the dynamics works at 
any energy. 

In [19L I20j . where a cylindrical lattice was considered, many points have been numerically or 
theoretically supported, for both homogeneous and disordered cases. For instance, starting from the 
fact that magnetic and kinetic energies are non-correlated observables, that the Boltzmann distribution 
is recovered at equilibrium and that the system results to be ergodic at all finite temperatures, it is possible 
to define the temperature by the averaged kinetic energy (E^): this is a link observable, and therefore 
a local quantity working as a temperature also in stationary states far from equilibrium, i.e. states forced 
by thermostats at different temperatures. 



In order to study the transient regime, we deal here with the same simple cylindrical geometry considered 
in [IHllSn]- It consists of a rectangular lattice of Lx x Ly sites, periodical in the Y direction and open 
in the X direction: the first and last columns, discrete circles of Ly sites, are therefore open borders, in 
contact with heat reservoirs constituted by supplementary columns (two for each side are enough). These 
reservoirs are kept at constant temperatures T\ and Ti by using the Metropolis dynamics at equilibrium 
with appropriate temperatures. In most cases, for our present exigences, only the thermostat at the left 
border is fixed, while the right border is open and virtually placed at infinity (i.e. very large Lx)- In the 
following we report and discuss the results obtained by means of Monte Carlo (MC) simulations performed 
on cylinders initially set at zero temperature, that is displaying an ordered magnetic configuration and 
zero kinetic energy on links. A given site i will be denoted by the couple {x, y}, where the integers x and 
y represent the distance from the reservoir and the position along the column respectively. Moreover, in 
order to focus on the main features of the transient, presently we limit to the homogeneous case, i.e. the 
coupling constants are set equal to 1 for any couple of neighboring sites. 

Time is measured in MC steps, each corresponding to a total extraction of Lx X Ly links, with 
consequent operations 1. and 2., as described in the previous section. In this way we get temporal series 
for the main observables, such as the average magnetization M(x,t) and the average energy E(x,t), for 
the x-th column at time t, being 



where o~ x ^ y is the value of the spin on site (x, y) at time t and (.) is the average over different thermal 
realizations (typically in our simulation we average over 10 4 thermal histories). Thus, as the simulation is 
onset, heat starts to propagate within the cylinder: the absolute average magnetization decreases and the 
average energies of columns close to the thermostat increase. 

The first part of simulations is devoted to the analysis of heat propagation, and the longitudinal size is 
taken rather large and equal to Lx = 300, so that heat can flow for long times without feeling the influence 
of the second border. In the following part we focus on finite-size effects and we consider relatively small 



3. Results 
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Figure 1. (Color on line) Main figure: Each curve, depicted in different color, represents the magnitude 
of the mean energy —E(x,t), as a function of the rescaled time t. Different columns of the system x have 
been taken into account. Inset: the energy E (red symbols) as a function of time t, compared with the 
power law ~ t 1 / 2 (black line). Numerical data can be fitted by a power-law with exponent 0.48; the small 
difference with respect to the expected value 0.5 can be ascribed to the existence of a long pre-asymptotic 
regime. All data represented in this plot refer to a system of size Ly = 300 initially ordered and subject 
to a thermostat at T = 10000. 



longitudinal sizes, that is Lx = 20. In any case, Ly is kept fixed and equal to 300; this large value is useful 
to avoid finite size effects and improves the statistics of the measures. 

3.1. Diffusion features 

In the presence of a static heat flux, the system behaves diffusively, as already noticed in [19 , 20 . The present 
simulations confirm this picture even in a transient regime. In particular, relevant physical quantities such 
as the total energy along the Y direction 

Lx 

E{t)=Y,E{x,t), 

x=l 

grows as i 1 / 2 , that is the characteristic law of diffusive processes; this asymptotic regime is reached only 
at very large times due to the extensive routines of the dynamics (see inset of Fig. I3.1[) . Moreover, 
Fig. 13. II evidences that focusing on the direction of heat propagation, the system can be studied within the 
framework of dynamical scaling which is typical of out-equilibrium systems. More precisely, the quantities 
(e.g. the local energy) depending on both time t and distance x can be plotted in terms of a single rescaled 
variable x — x/£(t), where i(t) ~ t 1 / 2 is the characteristic length of the process, so that data obtained 
at different spaces and times can be plotted together on a single curve. The fact that in our system £(t) 
grows as t 1 ! 2 confirms the diffusive behavior of the model. Of course, an equivalent scaling (which has 
been used in our plot) consists in introducing a rescaled time t = t/x 2 . Moreover, we remark that the 
scaling approach is very general and does not depend on the observable we have considered. Indeed we 
have checked that different observables, such as the magnetization and the correlation functions, display 
an analogous diffusive dynamical scaling. 

Finally, it is worth deepening the relation between the dynamical scaling introduced above and the 
generalized Fourier equation used in j!9j for the description of the model 

" E (3) 



dt dx \ dx _ 

where D(E) is a diffusivity factor which, in quasi-equilibrium states, is just the ratio between the 
conductivity and the specific heat, so that the standard Fourier equation is recovered. We recall that 
the previous equation is able to describe the heat propagation in very general conditions, independently of 
the existence of a local temperature [22]. Now, it is straightforward to verify that posing E(x,t) = E(x) 
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Figure 2. (Color on line) Log-log scale plot of the correlation function C n , for n = 1 (blue), n = 2 (red) 
and n = 3 (black), as a function of the rescaled time t = t/a: 2 . Measures have been performed on two 
columns, at distance 10 (•) and 20 (x) from the thermostat, respectively; due to time rescaling, data 
pertaining to different columns are collapsed. The green curves represent power-laws functions ~ t~ n / 2 
and provide a good fit for experimental data, over a proper time range. All data represented in this plot 
refer to a system of size Ly = 300 initially ordered and subject to a thermostat at temperature T = 10000. 



leads to 



xdE 
2~di 



d_ 

dx 



D(E) 



dE 
dx 



(4) 



i.e. the equation can be rewritten in terms of x only. Moreover, assuming that D{E) is a sufficiently 
smooth function writable as D{E) = Da + DiE, from Eq. [4] we can derive that the energy can be written 
as E(x) = Eq + Eix, namely E(t, x) — Eq ~ x/^/t (at least at long times). As we will see in the following 
subsection, this estimate is consistent with our numerical results. 



3.2. Correlations 

The transient regime presents interesting properties not only along the direction of heat flow but also along 
the orthogonal lines. In this framework, the simplest observable to consider is the two-point correlation 
function between spins of the same column at distance n, i.e.: 

1 Ly 

C n (x,t) = — y~]{<r x ,v(t)<rx,v+n(t)) ■ (5) 

Clearly, C\(x,t) is simply the negative of the energy E(x,t) of the x-th column, i.e. Ci(x,t) = —E(x,i). 
Since we are interested in the long-time asymptotic behavior, we consider the quantity 

C n (x,t) = C n (x,t)-C™, (6) 

where = (o- x ,yO~ x ,y+n) = (o-x,yO~ x +n,y) is the equilibrium two-point correlation function between spins 
at distance n in a 2-dimensional system evaluated at the same temperature of the thermostat. Therefore, 
if the system approaches to the equilibrium, we have lim t _ i . 00 C n (x, t) = 0, as it has been indeed verified in 
all our numerical simulations. 

In the Ising model outside the critical regime, when simply connected domains are prevalent, the 
correlation functions decay exponentially with the distance n, i.e. as exp(— n/X), where A is defined as the 
correlation length. At equilibrium, A depends on the temperature, while, in the coarsening processes, it 
grows as t 1 ^ 2 , due to the diffusion of the magnetic domains [TJ [3j [4] . 

Fig. 13.21 shows the behavior of C n (x,t) when the thermostat is at very large temperature, i.e. infinite 
for all practical purposes. In this case the equilibrium correlation function vanishes. The dynamical scaling 
picture is verified, since data evaluated at different columns (e.g. x = 10 and x = 20) collapse by using the 
rescaled variable t — t/x 2 , i.e. C n (x,t) — F n (t). Moreover, for large enough t, simulations evidence that 

F n {t) ~ r n ' 2 . (7) 

Therefore, the correlation function presents an exponential decay with the distance n and this is peculiar 
of the Ising model outside the critical regime. 

On the other hand, an unexpected and interesting feature is the slow decay of F n (t) as a function 
of t, consistently with the predictions drawn in the previous section for n = 1. While in the standard 
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Figure 3. (Color on line) Main figure: log-log scale plot of the correlation function C n , for n = 1 (blue), 
n = 2 (red) and n = 3 (black), as a function of the rescaled time t, for a system of size Ly = 300 
initially ordered and subject to a thermostat at temperature T = 5. Measures have been performed on 
two columns, at distance 10 (•) and 20 (x) from the thermostat, respectively; due to time rescaling, data 
pertaining to different columns are collapsed. The green curves represent power-laws functions ~ T ) 
and provide a good fit for experimental data, over a proper time range. Inset: Exponent a{n,T) obtained 
by fitting data from numerical simulations, as a function of n. Different temperatures for the thermostat 
are considered and compared, as explained by the legend. 



MC approach a slow equilibration is typical of freezing only (below the critical temperature) , here we can 
appreciate a similar behavior also during the warming phase. In fact, a rapid equilibration would involve 
an exponential decay in t. In particular, the heating process can not be considered concluded even for 
those sites whose distance from the thermostat is much smaller than £(t), i.e. the characteristic length of 
heat diffusion. From a different point of view, the very slow approach to equilibrium can be described in 
terms of correlation length A which, in our model, turns out to tend to zero logarithmically in i; indeed the 
correlation function can be written as: 

C n (x,t) = F n (t) ~ i- n ' 2 = exp 

We will show that the unexpected slow relaxation at high temperature emerges in open systems, while in 
the next section we verify that in small systems, where the boundary effects can not be discarded and heat 
can not flow freely, the thermalization process is exponentially fast at large enough times. 

The scaling function in Eq. [3] presents a very simple structure. Let us define g n (x) = F n (t) and expand 
g n (x) for small x (i.e. large times); we get 

oo 

g n (x) = J2A k , n x k . (9) 

fc=0 

Then, formula |(7J) can be recast as Ak, n = if k < n, giving a description, at infinite temperature, of 
the asymptotic behavior of the model in terms of the analytic properties of g n {x). Now, when the border 
is thermalized at a finite temperature, the main features characterizing the heat transport are preserved: 
as evidenced in Fig. 3, scaling still holds with a definite diffusive length, the correlation functions decay 
exponentially and the approach to equilibrium is slow since the scaling functions vanish as a power law. 
However, here the analytic form of F n (t) is not so simple; in fact we have F n ~ t~ a ( n ' 7 J, and fitting- 
procedures suggest that 

a{n 1 T) = a 1 {T)n + a Q (T) (10) 

as shown in the inset of Fig. 3. 

3.3. Dynamical regimes 



n log(i) 



= exp(— n/X). 



(8) 



The exponents characterizing the slow decays appear to be temperature dependent with a non-trivial 
behavior. However, Fig. 3 evidences that for high (practically infinite) temperature, a(n,T) — n/2, while, 
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for low enough temperatures, a(n, T) ~ 1/2 independently of n. This suggests that two simple (asymptotic) 
behaviors are present in the model: for infinite temperature the scaling function should be given by 

hence, at infinite temperature, the dynamics is characterized by indefinitely shortening correlation length 
A(f) = 2/log(i). This is clearly consistent with the fact that at infinite temperature there is not an 
equilibrium characteristic length, since the correlation function between any couple of sites is vanishing. 
On the other hand, at small temperature, F n (i) decays as i~ x / 2 . However, since the exponential decay of 
C n (t) as a function of the distance n is an expected feature, we have 

F n {t)~t- l ' 2 eM-n/\') (12) 

where a characteristic length A' naturally emerges. We remark that at finite temperature the equilibrium 
correlation functions decay with a characteristic length X eq (T). Since a similar behavior is expected to 
be reached asymptotically also out of equilibrium, we expect that, eventually, A' ~ X eq (T). Therefore the 
above regimes (|11I12I) provide a reasonable description for infinite and finite temperature respectively. 

In this perspective, the complex behavior characterizing intermediate temperatures would result from 
an interplay between these regimes. In particular, when the thermostat is at a finite temperature T, 
F n (t) evolves according to Eq. (fTTj) until the correlation length A(t) becomes comparable to A' ~ X eq (T); 
afterwards an evolution described by (|12l) takes place. At very large temperatures A' ~ X eq (T) is almost 
zero and, in the considered time regimes, X(t) 3> A', therefore only the infinite temperature evolution (llip 
is observed. At low temperatures the correlation A' is so large that A(f) becomes quickly comparable to 
A' and F n (t) evolves according to (|12[) . Finally, the temperature dependent exponents a(n,T) observed at 
intermediate temperatures should be the result of the crossover between the two regimes which is ongoing 
in the considered time window. 



F„(t)~i-"/ 2 = exp 



3.4- Finite-size effects 

Let us now discuss the finite-size effects in the processes of heat transfer from the thermostat to the 
system. Finite-size effects can clarify some features outlined in the previous analysis, in particular the 
mutual influence between the transport regime and the way the asymptotic values are approached. 

Fig. 4 shows the behavior of Ci(x, t) as a function of time for the column x = 10 and a size Lx = 20. 
Differently from the behavior depicted in Fig. 13.11 here the dynamics speeds up, approaching quickly the 
equilibrium value, for large enough times, i.e. when £(t) is comparable to the lattice size. Therefore, the 
presence of very long transients described by a power law seems to be typical of infinite systems, while 
finite systems thermalize in a rapid (possibly exponential) way. This differences ensure that the dynamics 
may be used for the study of the stationary properties of finite systems as it has been done in [T21 [201 HI] ■ 
Moreover, Fig. 4 evidences that the characteristic time for equilibration increases with the temperature 
becoming very large at T = oo. 

4. Conclusions 

The dynamical insight on the warming process in the Ising model we have studied has two main aspects: 
first, a direct confirmation of such expected features as the diffusive behavior in the heat propagation, or 
the exponential decay of correlations in the transversal direction; second, a clear indication about the role of 
heat flux in keeping the time behavior of correlations well distinct from their time-asymptotic "equilibrium" 
value. Indeed, in the indefinite propagation taking place when the longitudinal size is much longer than the 
observation window, the approach to the limit obeys a power law in time. Remarkably, a slow relaxation 
seems to be present for all temperatures. In fact, one could in principle conjecture the existence of a special 
temperature (possibly related to the critical one) separating a slow and a fast dynamical regime; yet the 
persistent slowness of the dynamics proves the non-existence of such a dynamical threshold. In particular, 
the explanation of the non-trivial behavior of the temperature dependent exponents, in terms of finite and 
infinite temperature regimes, confirms the fact that other intermediate thresholds do not exist. 
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Figure 4. (Color on line) Correlation function Ci as a function of time t, for a system of size Ly = 300, 
Lx = 20 and x = 10, subject to a thermostat set at different temperature, each temperature being 
represented by a different color as explained by the legend. 
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